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Abstract. We examine the quahtative features of parton production through 
materiahzation in heavy-ion coUisions within perturbative QCD, and estimate 
the magnitude of the resulting parton density created during the early stage 
of the coUisions. The implications for "anomalous" J/ip suppression observed 
in Pb-I-Pb collisions at the CERN SPS are discussed. We argue that the A- 
dependence of absorption of J/ip by (partonic) comovers is steeper than as- 
sumed in most phenomenological models, because the absorption process is 
dominated by quasi-perturbative QCD interactions. Our argument is sup- 
ported by results recently obtained in the framework of the parton cascade 
model. We predict significant "anomalous" suppression for Pb-|-Pb collisions 
at the CERN-SPS, but not for S+U coUisions. 
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1. Introduction 

The creation of high-density QCD matter consisting of quarks and gluons that are 
being liberated in ultra-relativistic heavy-ion collisions, is one of the most funda- 
mental issues for investigating the new physics regime of deconfined parton matter 
and the quark-gluon plasma. In this context, the question whether the NA50 data 
[01 on J/tp production in Pb + Pb collisions at the CERN-SPS signal the presence 
of a novel nuclear suppression mechanism is a matter of heated controversy. Some 
theoretical papers |^-^ have argued that the additional suppression observed in Pb 
+ Pb, compared with S -I- U, is "anomalous" , i.e. provides evidence of a new mech- 
anism, most likely dissociation of the cc quark pair in a dense deconfined medium, 
or quark-gluon plasma 1^. Others [^-^ have argued that "conventional" mecha- 
nisms, such as absorption of J/ip on light hadrons emitted at the same rapidity, can 



0231-4428/97/ $ 5.00 
©1997 Akademiai Kiado, Budapest 



2 



K. Geiger and B. Miiller 



explain the increased suppression in Pb + Pb collisions. 

It is not our goal here to add another phenomenological analysis of the NA50 
data to the existing ones. Rather, we will make an attempt to analyze on the 
basis of the established concepts of perturbative QCD, whether Pb + Pb collisions 
{Pbeam = 158 ^ GeV) are indeed so different from S + U collisions at almost the 
same energy {Pbeam — 200 A GeV) that the former system might result in the 
formation of a thermalized, color-deconfined state whereas the latter does not. Our 
reasoning is based partly on qualitative arguments first annunciated a decade ago 
and partly on quantitative results recently obtained in the framework of 
the parton cascade model of high-energy nuclear collisions |l2,0|. 



2. A-dependence of parton production 

We begin with the qualitative arguments. Energy deposition at large transverse 
momentum p± in hadronic collisions can be described in the framework of pertur- 
bative QCD through jet or mini-jet production. The crucial question is how low in 
p± this description remains reliable. There is an increasing amount of evidence fl^ 
that the perturbative domain of QCD extends down to (spacelike) momenta of the 
order of 1 GeV/c. Contrary to hadron collisions at presently accessible energies, 
in heavy-ion collisions the applicability of perturbative methods involves an addi- 
tional scale namely the parton density in the colliding nuclei. For large nuclei even 
at moderately small values of the Bjorken variable x, the Lorentz contraction of the 
nucleons may lead to overlapping partons with the net effect that parton densities 
appear to be suppressed as compared to isolated hadrons. Hence, not only the 
momentum transfer p± of mini-jet production, but also the nuclear mass number 
A determine the perturbative regime of parton dynamics. 

In order to estimate the A-dependence of the density of partons produced in 
heavy nucleus collisions and to pin down the point of critical density where the 
parton densities begin to deviate, we follow ref. ]lO| by considering the number of 
mini-jets produced in a nucleus- nucleus collision A + B at impact parameter b [ p5[ : 

""^^2^ " 2rAB(6) 2^ J,^—X fa/A{x ,Px) Xfb/B{x,p^) ^ (1) 
ab -L 

where TAsib) = / d'^sTA{s)TB{b — s) with the nuclear profile functions Ta and 
Tb normalized to unity. The functions fa/ a and fb/B are the distribution of par- 
tons a,b = g,qi,qi in the colHding nuclei A and B, which depend on the energy 
fractions x, x' and the transverse momentum p± of the partonic scattering with 
cross-section aabix, x' ,p\). For purpose of simplicity, we consider now symmetric 
collision systems, i.e. A — B, but the arguments below generalize straightforwardly 
to asymmetric systems A ^ B. 

To proceed, we recall that for sufficiently large values of p± and not too small 
values of x, the nuclear parton distribution functions fa/A are approximately given 
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by 

fa/Aix,?^) ~ Afa/N{x,pi), (2) 

where fa/N is the number density of parton species a in the nucleon N. However, 
as X and p± are decreased, the parton density in a heavy nucleus is increasingly 
reduced ("shadowed") because partons from different nucleons overlap spatially 
and begin to fuse. Since the small- a; regime is strongly dominated by the rise of the 
gluon density with quarks being negligible, the behavior of the gluon distribution 
fg/N{x,p\) with decreasing x and may be viewed as representative indicator for 
shadowing effects to become significant. For gluons, the characteristic momentum 
Pent where this happens is determined by the number of gluons per unit area in the 
transverse plane of the colliding nuclei at this momentum scale pO[ | , 

2 9 Axfg/N{x,pl^.^) 
Pcrit = Ca ^1 , (3) 

where Ca = 3 is the value of the quadratic Casimir operator of SU(3) in the adjoint 
representation, and Ra = 1.2 A^^^ fm is the nuclear radius. Notice that p^^^^ enters 
as argument of fg/pj on the right side; therefore is a non-trivial implicit equation. 
Nevertheless, we may get a rough estimate of the condition (^) by taking for instance 
Pcrit = 1 GeV and xfg/N{x, IGeV^) = 3 at a; « 0.1 with as ~ 0.3 at p^rit' which 
case one finds that shadowing corrections would become significant only for nuclei 
with ^ > 240. 

Since Ra ^ A^^^, the critical transverse momentum grows approximately as 
Pcrit ~ A^^^. The effect of shadowing is to reduce the number of mini-jets with 
transverse momenta p± < Pcrit- For larger values of p_L, the mini-jet yield grows as 
A'*/^, because fg/A ~ ^fg/N and TAAib) ^ R'a' A^'^^'^. However, approximating 
the glue-glue cross section by its small angle limit 

dagg_ _ {asCAf-^ , .s 

dpi ^ 2p\ ' 

one finds that the integrated mini-jet yield for all momenta p^ > pcrit only grows 
linearly with A: 

This agrees with the result of Blaizot and Mueller who found that for two 
heavy nuclei the total yield of gluon mini-jets per unit of rapidity is 

dN{h) TAAjb) 

The same considerations apply when the mini-jet production is mostly due to quark- 
quark or quark-gluon scattering, because the differential parton cross sections all 
have the same functional form (^) at small angles. 
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To summarize, we have the following picture describing the mini-jet production 
at central rapidity in nucleus- nucleus collisions: Above a critical transverse momen- 
tum Pcrit, which itself grows as A^^^ , the mini-jet multiplicity grows like A'^/^. The 
total mini-jet yield, as well as the multiplicity in the saturation regime below Pcrit, 
only grows linearly with A. This result implies that the area density of mini-jets 
grows as A^/'^ above Pcrit, but only as A^^^ below Pcrit- 



3. Implications for J/ip suppression 

We now return to the problem of the observed strong suppression of Jf-ip production 
in Pb -|- Pb collisions. In particular, we want to address the question of the A- 
dependence of the absorption of the J/V' by comoving produced hadronic matter. 
In the framework developed by Bhanot and Peskin [|l6[, interactions between light 
hadrons and deeply bound, heavy quarkonium status, such as the J/ip, are mediated 
by short-range color dipole interactions. The relevant momentum transfers for 
these interactions lie in the perturbative domain of QCD, p]_> 1 GeV'^. 

The key element in the following arguments is the fact that, for most nuclear 
collision systems at the relatively low center-of-mass energy of the CERN-SPS ex- 
periments, perturbative mini-jet production is largely due to quark-quark scattering, 
because the typical value of Bjorken-x probed is (x) ^0.1, where the gluon density 
is small. This implies that shadowing effects are negligible, i.e. the relation (^) 
holds. Furthermore, according to the critical momentum Pait lies below the 
perturbatively accessible range of momenta pj_ > I GeV at the SPS energy, and 
may just begin to reach into it for the Pb -I- Pb system. If this is correct, then all 
mini-jet production involving momenta > 1 GeV lies safely above Pcritj where eq. 

holds. As immediate consequence, the density of comovers which can effectively 
interact and absorb the J/tp at the SPS grows like A^/^, or (^1^2)"'^^'^ in asym- 
metric collisions. This is a much stronger A-dependence than naively expected and 
embodied in most comover suppression models It also implies a stronger im- 
pact parameter or Et -dependence of comover suppression than predicted by existing 
models. 

Quantitative support for our speculations comes from recent calculations pT| ] 
of secondary particle production in the parton cascade model This model is 
based on perturbative parton scattering with dynamic, medium dependent cut-off 
mechanisms of the infrared divergences of QCD. The model also incorporates a 
density dependent clustering mechanism for the parton-to-hadron transition. The 
calculations predict that the energy density e produced by scattering partons (mini- 
jets) at central rapidity grows by a factor of more than 2 between S -I- U and Pb -I- 
Pb, compatible with the scaling law e ~ {AiA2)^^^ , but much faster than expected 
in the usual Glauber model approach which predicts e (^1^"^ + ^2^"^); in which 
case e increases only by 30 %. 

In Pb -|- Pb collisions the initial partonic energy density ep is predicted to reach 
5 GeV/fm^ at times r < 1 fm/c in the comoving reference frame (see Fig. 1). For 
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a thermalized gas of free gluons and three flavors of light quarks this corresponds 
to an initial temperature Ti « 230 MeV, clearly above the critical temperature 
Tc ~ 150 MeV predicted by lattice-QCD calculations shown in Fig. 2. On the 
other hand, the parton density ep w 2 GeV/fm'^ predicted for S + U collisions (Fig. 
1), just lies at the upper end of the transition region in the equation of state from 
lattice-QCD. 

The temperature Tb required for dissociation of the J /ip bound state due to 
color screening is known to be higher than namely To ~ 1.2 = 180 — 200 
MeV. The parton cascade model results are therefore compatible with the experi- 
mental finding that there appears to exist no significant comover-induced suppres- 
sion of J /ip in S -|- U collisions, but a large and strongly impact parameter dependent 
effect is observed in Pb + Pb collisions. 

In conclusion, we have argued that the density of comovers relevant for J /ip 
suppression in nucleus-nucleus collisions grows more rapidly as function of nuclear 
mass than assumed in previous studies based on the Glauber model. Recent results 
obtained in the framework of the parton cascade model are therefore consistent 
with the interpretation that the "anomalous" J/ip suppression observed in Pb -t- 
Pb collisions is caused by color screening in a dense partonic medium. The same 
calculations predict that there is no such effect in S -f U collisions. It must be 
stressed that these model results do not contradict the experimental data collected 
at the SPS, as one might suspect, since it is usually claimed that those are fully 
consistent with the Glauber picture. In fact, the detailed study of particle spectra 
in Ref. |ll|, comparing model and experiment, shows that the data are compatible 
with a parton cascade picture. In particular, the magnitude of transverse energy 
production visible in the final-state hadron spectra, is in fair agreement, although 
the maximum achieved energy density during the very early stage in Fig. 1 is 
substantially larger than previously estimated. 
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Fig. 1. Time evolution of the energy density e of the partonic matter in the central 
slice of the collision systems S + U and Pb + Pb at CERN-SPS beam energies of 
200 A-GeV and 158 A-GeV, respectively. The model calculations were obtained 
with the MC implementation |17| of the parton cascade model [^^. 
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Fig. 2. Equation of state for two-flavor QCD from lattice calculations by the 
MILC collaboration |Q , showing energy density e and pressure p as a function of 
temperature T. Indicated in the figure are also the energy densities corresponding 
to several values of the temperature T, as well as the location of the " melting" points 
of the J/ip and T states, where these bound states disappear due to dissociation 
0- 



